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0.1 INTRODUCTION 



The beautiful experiment measuring the muon anomalous magnetic moment (g-2) now in progress 
at Brookhaven (l| has already attained a precision that challenges our ability to calculate the 
expected result of the experiment. The uncertainty of the Brookhaven g-2 measurement is now 
about 1.5 in units of 10 10 a /i ( where = ^-5-), with 1 year of data yet to be analyzed according 
to reference 1. Weak and electromagnetic contributions to g-2 can be calculated according to a 
well- understood theory with about 1/30'th of the uncertainty of the measurement^], as is noted 
in the reference. The difference between that calculation and the measurement is 71.8 with the 
uncertainty of the experimental result. The difference is ascribed to hadronic effects (see |3| for a 
recent review). 

The significance of the new result is, in my view, that it measures the hadronic contribution 
to the muon g-2 with a precision of the order of 1 %. The precision may increase by about an 
order of magnitude when the experiment is completed. There is no theory of hadrons that predicts 
hadronic interactions so precisely. The appropriate theory should presumably be grounded in 
quantum chromodynamics (QCD), which does not yet lend itself to such precise predictions. There 
are phenomenological arguments which relate the hadronic contribution to g-2 to other experiments 
and provide the basis for current estimates which range § from 69.2 to 72.5, each with a quoted 
uncertainty of the order of a percent. 

The quoted uncertainties of current estimates reflect only the uncertainties of the experimental 
data used for phenomenological estimates and the computational uncertainties of the estimates. 
They do not, in my view, adequately reflect the uncertainties in the theories underlying the phe- 
nomenological estimates, as I shall explain in the remainder of this note. 



0.2 DESCRIPTION OF THE HADRONIC CONTRIBUTION 

In the notation of Sommerfield the muon propagator G mu relates to a mass function M as 

G^ = 7 n + M (1) 

where = (p — eA)u, p is the muon momentum and A is the potential for a weak constant external 
field. The anomalous g-factor occurs as a term proportional to (XnvFav in the mass function M. 
That function is calculated from the equation 

M = m + ie 2 1 G mu TQ (2) 

involving the photon Green's function Q and the "dressed photon-muon vertex " T. 

The photon Green's function is assumed ||] to obey a Kallen-Lehman representation Q. Ac- 
cording to that representation each element Qn U may be written as a dispersion integral 

^) = ^-[^-^]f (3) 
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The weight function s(x) is proportional to the sum over possible intermediate states produced by 
the electromagnetic current, according to the equation 

k A s(-k 2 ) V = \ ]T < I j M (0) I a >< a I j„(0) | > 5\a - k) (4) 

where j is the electromagnetic current operator. The hadronic contribution to Q is then given by 
that part of s(x) for which the states | a > consist of hadrons. The matrix elements < | j „ \ a. > 
are the amplitudes for production of the state | a > from the vacuum by the electromagnetic 
current. The corresponding contributions to s(x) are therefore proportional, to lowest order in 
the fine-structure constant, to the cross-sections for producing states a in e + e _ collisions. It 
would accordingly seem that such measured cross-sections could be used to calculate the hadronic 
contributions to Q. 

The components of the vertex function T contain terms corresponding to the production of 
a virtual hadron that couples to the muon. The lowest order such terms (in powers of the charge 
e) - the so-called "light-light scattering" terms - are of order e 6 . Power counting suggests that 
such terms contribute to the anomalous muon moment at a level of about ^jth of the uncertainty 
in the Brookhaven measurement. There is, however, no theory comparable to electroweak theory 
for calculating these terms, as Hayakawa and Kinoshita have emphasized Q in connection with 
their own estimates. Those estimates are not severely tested at the present level of experimental 
accuracy. 



0.3 UNCERTAINTIES IN HADRONIC-CONTRIBUTION ES- 
TIMATES 



The Kallen-Lehman representation [|7]] for the photon propagator is of the nature of a dispersion 
relation. The dispersion relation is a statement that Qp L)U {k 2 ) is a real analytic function of the 
variable k 2 , cut along the negative (that is, timelike) k 2 axis Q. It is also, however, a statement 
that Q is a polynomially-bounded function for large, complex k 2 . I contend that the latter statement 



is unprovable [10] [|ilf| . There does not, in any event, appear to be any existing proof that Q obeys 
such a boundedness property. Any attempt to assign a numerical uncertainty to a calculation 
based upon a Kallen-Lehman representation for the photon propagator must therefore be regarded 
as speculative. 

There are recent estimates of the hadronic contributions to the muon g-2 that are within a 
few per cent of values that one would deduce from the new experimental result Q. Is this near 
agreement evidence for the polynomial boundedness of G? I present in the next section models, 
by way of counter-examples, where such near agreement can occur even with functions that grow 
exponentially. Polynomial boundedness cannot, therefore, be deduced from the present data. Un- 
certainty estimates based upon dispersion-relation calculations of Q must be understood as lower 
limits because the theories underlying the calculations are incompletely defined. 
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0.4 Dispersion Relations for Functions Having Exponential Growth 



I present two models to show that an attempt to calculate a function by a dispersion relation, when 
the assumption of polynomial boundedness is incorrect, need not lead to a result that is wildly 
different from the correct result. The fact that dispersive calculations of the hadronic component 
of the muon g-2 are close to the values deduced from experiment is, therefore, not evidence that 
such calculations can give correct results. 

Both of the models happen to involve functions related to Bessel functions. The essence of the 
demonstration, however, may be made with the two elementary functions natural logarithm and 
exponential. Consider, for this purpose the following function Q(z) of a complex variable z: 

Q{z) = ln(-z) + e~ bz (5) 

where b is some real positive constant. Q{z) is an analytic function of the complex variable z except 
for logarithmic branch points at and oo and an essential singularity (of the exponential) at oo. 
Q is made single-valued by a cut along the positive real axis, giving Q an imaginary part equal to 
— iir as the real axis is approached from above. 

The point of the models is that they involve functions that, like Q, have an essential singularity 
and are therefore not polynomially bounded. One may nevertheless try write a dispersion relation 
for Q by integrating around the cut, because the imaginary part, being a constant in this case, 
is well-behaved. The result of the dispersion integral will be to recover only the In part of Q. 
The remaining part of Q, the exponential, will be negligible (compared with the logarithm) almost 
everywhere along the positive real axis, depending upon the size of the parameter b. The erroneous 
assumption that Q satisfies a dispersion relation therefore makes only a small, model-dependent 
error almost everywhere along the positive real axis. 



0.4.1 A Hankel- Function Model 

In the first model I assume that the correct function that I am trying to calculate is the Hankel 
function[|l2| H^\z). The imaginary part of the Hankel function is the Neumann function Yq{z). 



Both functions emulate a coulomb-like singularity at z = 0, albeit only logarithmically, so a numer- 
ical integration of a dispersion relation must be cut off at small z. The function Hq \z) behaves 
for large, complex z like ^7=, and is not, therefore, polynomially bounded. 

The model is unrealistic in the sense that the imaginary part of the "true function" is oscilla- 
tory, unlike the strictly positive cross-section data that is used as input in the photon-propagator 
calculations. That particular unrealistic feature is corrected in the second model. 

The game, then, is to pretend that we know from experiment the imaginary part of the "correct" 
function. We then put the "experimental" data into an unsubtracted dispersion relation under the 
erroneous assumption that the "correct" function vanishes sufficiently rapidly at large \z\, cutting 
off the integration at some small value of z. We want to determine how much our answer deviates 
from the correct answer which is, in this case, the Bessel function Jq{z). 

The answer can be determined analytically, because 
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1 f°° Y (x)dx 2 

- Jo{z) , 0-1,0(2) (6) 



TT JO X — Z 7T^ 



where the principal value of the integral is intended and S—io(z) is a Lommel function whose 
explicit form has been given by Watson [|13||. It is a solution of the inhomogeneous Bessel Equation 



S" + -S' + S = \ (7) 



z 2 



subject to the condition that it behaves asymptotically like \. It grows for small z like ln 2 (z), 
but rapidly becomes small as z increases. Fig. 1 compares the right hand side of Eq. ||, shown by 
the solid line, with the "correct" answer Jo(z), shown by the dotted line. It is apparent that the 
error from using misusing the dispersion integral can be made quite small if small values of z are 
excluded, as they are in "real life" when coulomb and infra-red effects play a role. 

0.4.2 A Strictly Positive Imaginary Part 

A more realistic model would require that the simulated experimental data be given by a function 
that is strictly positive, like a cross-section. A convenient choice is the modified Hankel function 
Kq(z) which approaches zero asymptotically for large z like ^-j= and grows logarithmically for small 
z like —ln(z). For z approaching the positive real axis from above, there is a function f(z) given 
by 

l r°° 11 

f(z) = - K (x)dx( —) = ir(L o (z)-I (z)) + tK (z) (8) 

TT Jo X — Z X + Z 

Here Lo is the modified Struve function and Io is the modified Bessel function, each of order zero. 

The complex function f(z) does not satisfy a dispersion relation. This is because the integral of 
the function around a closed loop at large z cannot be made arbitrarily small, since its imaginary 
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FIG 2 



part, -Zt"o (^) 5 grows without bound when | < arg(z) < 4p The imaginary part of /(z) again 
represents the simulated experimental data. 

The game, as before, is to put the simulated experimental data into an unsubtracted dispersion 
relation and compare with the correct answer, given in Eq |8[ The result is shown in Fig. 2 which 
displays the real part of the dispersion result as a solid line and the real part of /(s) as a dotted 
line (the imaginary parts are, of course, identical). It is evident that for values of z greater than 
about unity, the dispersion result is not greatly different from the "correct" value. 

0.5 Conclusion 

Hadronic contributions to g-2 originate, according to our present understanding, from electromag- 
netic excitation of virtual quark pairs from the vacuum. The task of calculating the effects of such 
excitations therefore falls within the realm of Quantum Chromodynamics ("QCD"). Such calcula- 
tions have not yet been done. It has been tempting to suppose, as an alternative, that dispersion 
relations using empirical data could be used to predict the results of QCD calculations, at least for 
the vacuum polarization corrections to the photon propagator. 

The difficulty with the supposition is that a dispersion relation involves an assumption about 
the behavior of a function for large values of its argument, namely, that the function is at least 
polynomially bounded. In the case of the photon propagator there is no way to justify that as- 
sumption. 

It might be tempting to try to evade this lack of justification by reversing the burden of proof, 
asking: "What is the physical origin of the analogous terms in the physically relevant case? In 
particular, if Exp[-z] — > Exp[-E/E0], what is the meaning and value of the energy scale EO?" 

Because the presence of an essential singularity at infinite energy, as I have modeled with an 
exponential term, is purely conjectural, any mass, such as the pion mass, the QCD scale or a 
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spontaneously generated mass [14] is a possible candidate for the energy scale EO. One doesn't see 
such terms in the usual discussions of analytic properties of amplitudes because those discussions 
are customarily guided by perturbation theory [O] where such terms cannot occur. But we 



know of models where perturbation theory breaks down. [14]. And essential singularities are a 
common feature of solutions of differential equations such as the relativistic Schroedinger equation, 
cf. 0. 

Proof of the existence of an essential singularity would, of course, end the discussion; dispersion 
relations would not provide an algorithm for calculating the hadronic contribution to vacuum 
polarization. That is not the intent of this note, however. The intent is merely to point out that 
there is, as yet, no mathematical justification for the use of dispersion relations for that purpose, 
and the accuracy of the resulting predictions must therefore remain uncertain. 

What I have shown here is that violation of the assumption of polynomial boundedness can, 
in some circumstances, lead to results that are close to the correct results. There is, however, 
no apparent a priori way to estimate the size of the error resulting from the violation. There is 
consequently no way to estimate the error involved in using a dispersion relation to calculate the 
hadronic contribution in question. [17] 

The Brown, et a/.[|J measurement of the muon g-2 is therefore significant because it measures 
the hadronic contribution to g-2 with unprecedented precision. That measurement stands as a 
challenge to our understanding of the quark structure of the vacuum. [17] 

I am indebted to Geoffrey Bodwin, Tim Tait, Carlos Wagner, Alan White, and Cosmas Zachos 
for some helpful discussions, to Professor T. T. Wu for his assistance in identifying reference Jl4| , 
and to the (anonymous) referee, whose question led me to become aware of references ] pj] , [14|, [|16| 
and several related papers. This work was supported by the U.S. Department of Energy, Division 
of High Energy Physics, Contract W-31-109-ENG-38. 
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